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ABSTRACT 

Aims. We explore planet formation in binary systems around the central star where the protoplanetary disk plane is highly inclined 
with respect to the companion star orbit. This might be the most frequent scenario for binary separations larger than 40 AU, according 
to Hale (1994). We focus on planetesimal accretion and compute average impact velocities in the habitable region and up to 6 AU 
from the primary. 

Methods. Planetesimal trajectories are computed within the frame of the restricted 3-body problem determined by the central star, the 
companion star and massless planetesimals. Relative velocities are computed and interpreted in terms of accreting or eroding impacts. 
Results. We first show that, for binary inclinations higher than 10 degrees, planetesimals evolve, at a first approximation, in a gas-free 
environment. Planetesimal accretion is confined around the central star to a region determined by two main parameters, firstly by the 
mutual inclination between the binary plane and the disk, and, secondly, by the binary eccentricity. 

Conclusions. The onset of large mutual inclinations between planetesimals due to the nodal randomization causes an increase in 
the relative velocity. The chances for a successful planet accumulation process depend on the balance between the timescale for 
node randomization and that of planetesimal accretion. When the binary semimajor axis is larger than 70 AU, planet formation 
appears possible even for eccentric binaries (up to 0.4). For lower binary separations the region where planetesimals accumulate into 
protoplanets shrinks consistently. When the mutual inclination between the binary plane and that of the planetesimal disk is larger 
than 40° the Kozai mechanism strongly inhibits planetesimal accumulation. 

Key words. Planetary systems: formation; Celestial mechanics; Methods: numerical 

1 . Introduction domly inclined respect to the binary planet. As a consequence, 

the inclination between the binary's orbital plane and the cir- 

Planetary formation in binary systems is a complex issue, cumprimary disk is a parameter that has to be taken into account 

since each step of the process can be affected in different when stu d y ing planetesimal accretion, at least for systems with 

wa ys by the companion per t urbations. Recent numerical stud- a > 30 - 40 AU 

fes (Marzari & ScMiJ^^p^LMM M W e will focus in this paper on binaries with intermediate sep- 
£ ; ^S&dekooper^tSpQ08; Xie &Zhm4£QQ8D have shown fa ^ ^ ^ ^ indination * 
£ . that one stage is particularly sensitive to the presence of the sec- a free meter similarly to the studies for the coplanar case, 
ondary star: the initial accretion o f kilometer-siz e^anetesimals ±e main QUtcome wg ^ interested on is the ; t velocit dis . 
(a review on this topic is given in| Haghighipouij (|2009|)). Indeed, tribution within the planetesimal population, since this parame- 
the coupled effect of secular perturbations of the companion star ^ CQmoh ^ fete of lanetesimal co ii isions , either accretion 
and friction due to gas in the nebulae induces a size-dependent Qr eroskm For a si ificant mutual inc i ination betW een the bi- 
phasing of orbits which may lead to high impact velocities. This Qrbital , ane and ^ digk of planetesimals embedded in the 
could slow down or even halt the accretion process even in the us ^ ^ forced inclination due to the companion star 
terrestrial planet region for a wide range of binary separations, mi ht bg much mQre effectiye ^ increasing the relative veloc- 
i.e., up to a, ~ 50 AU forhi^heccent^ systems (see for ex- and haki ^ formation As the pi anet esimals decou- 
ample Figs.8 and 9 in| Thebault et all |2006|). ple from the gaseous difjk and evolye grav i tat i ona n y , they wou i d 
However, these studies are based on the assumption that the feel the binary perturbation and move into inclined orbits. The 
planetesimal disk is coplanar to the stellar orbit. Even if this as- perturbations of the companion star leads to a progressive ran- 
sumption might intuitively seems reasonable, a systematic study domization of planetesimal node longitudes, starting from the 
by|HaleJ dl994j) on binary systems with solar-type components outer region of the disk where the secular periods are shorter, 
suggests that the spin of the two stars is aligned only for binary xhe planetesimal disk gradually evolves into a cloud with an an- 
systems of 30-40 AU or less. Beyond this distance, the primary's gular ope ning equal to twice the mutual initial inclination of the 
equator, and thus a putative planetesimal disk, appears to be ran- disk reS p ec t to the binary plane. We explore in this paper the 
effects of the nodal randomization on the mutual relative veloci- 
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We also estimate the minimum inclination below which planet 
formation may occur in spite of the binary inclination. 

The paper is organized as follows: in Section 2. we show 
that the planetesimal dynamics perturbed by the companion star 
keeps the swarm out of the gaseous disk for most of the orbital 
period. This makes gas drag a negligible perturbation. In Section 
3. we describe the numerical model used to compute the plan- 
etesimals relative velocities. Section 4. is devoted to the anal- 
ysis of the impact velocities for different binary parameters. In 
Section 5. we derive limiting inclinations for accretion at differ- 
ent binary separations. Finally, in Section 6. we summarize our 
results. 




2. Decoupling between gaseous disk and 
planetesimals 

Most rece nt studies of planetesimal accretion in a binary envi- 
ronm ent dMarzari & Scholl 120001: iThebaultet all I2006L 120081 
2009; Paardekoope r et al.L 120081) have focused on the influence 
of the gaseous component of the disk on particle dynamics. 
However, the implicit assumption that the planetesimal swarm 
is embedded in the gas disk is only valid if the disk is copla- 
nar to the binary orbital plane. In this case, planetesimals feel a 
steady gas drag and have their orbital evolution significantly af- 
fected by frictional forces. However, if the companion star is on 
an inclined orbit with respect to the disk mid-plane, the situation 
is dynamically more complicated. Three possible scenarios can 
be envisaged for the interactions between gas and planetesimals: 

- Planetesimals form within the gas disk which remains a 
long-lived coherent entity in spite of the binary perturba- 
tions. Numerical simulations with constant viscosity and a 
polytr opic equation of state performed by lLarwood et alj 
(1996) with an SPH code suggest that a disk perturbed by 
an inclined companion star maintains a coherent structure if 
the Mach number is lower than 30. It behaves like a rigid 
body preceding at a rate a> p given by: 



15M S R 3 

dir. -cos(i m )Q,(R) 

' 32M p D 3 



(1) 



where M s and M p are the masses of the secondary and pri- 
mary star, respectively, D is the radius of the circular orbit of 
the binary, i m is the mutual inclination between the disk and 
the binary orbit, Q(R) is the keplerian frequency and Rj is 
the disk radius. This equation is derived under the simplified 
assumption that the disk has a constant density, but it is in 
general a good approximation to more general cases. In this 
scenario, when the planetesimals reach the size (1-10 km in 
diameter) for which they evolve under the dominating grav- 
itational force of the two stars, they leave the disk plane be- 
cause of the forced component in the inclination. Their orbits 
move in the binary orbital plane and their nodes circulate at 
different rates, depending on their semimajor axis. Gas drag 
is probably not a significant perturbation in this scenario, 
since it affects planetesimal evolution only in the fraction of 
time during which they cross the gaseous disk plane. This is 
clearly illustrated in FigQ] where we show the projection of 
the planetesimal positions with respect to the gaseous disk 
when the inclination of the binary orbital plane is i m = 20° 
with respect to the initial disk plane. Planetesimals spend 
most of their time out of disk where the gas density is neg- 
ligible. According to our simulations, for i m = 30° planetes- 
imals spend on average only 9% of their time within one 



X(AU) 

Fig. 1. Planetesimal positions (red squares) in the x-z carte- 
sian plane after 1 x 10 5 yrs from the beginning of their gravity 
dominated evolution. The gaseous disk (assuming a scale height 
h = 0.05(r/At/) (5 / 4) AU) is shown by green dashed lines. The 
orbital plane of the binary is in the x-y plane and the disk is 
inclined of i m - 20°. 



scale height of the gaseous disk. This fraction increases to 
13% when i m = 20° and to 27% when i m = 10°. As a conse- 
quence, we estimate that for i m > 20° gas friction can be, to 
a first approximation, completely neglected when computing 
planetesimals orbital evolution, while the i m = 10° case ap- 
pears as a limiting value below which gas friction has to be 
taken into account. 

- The gaseous disk begins to warp and it is disrupted by the 
binary perturbations. It loses c oherence and the gas is dis- 
persed in space. According to lLarwood et al.l (I 19961) such 
disruption by differential precession might affect extremely 
thin disks. Also in this case, the planetesimals would evolve 
in a gas-free environment. If the disk disrupts before the 
planetesimals detach from the disk then this would be the 
most significant gas free case, where planet formation would 
start from a disk of solid material made of small planetesimal 
precursors which would evolve under gravity only. 

- As in the first case, the disk remains coherent but it relaxes 
to the bina ry plane on a timescal e comparable to the viscous 
timescale (lLarwood et al.L 119961) . If the process is fast due 
to a high viscosity of the disk, kilometer-sized planetesimals 
have not the time to form before the disk relaxes to the bi- 
nary plane. Planetesimals would then grow when their orbital 
plane, and that of the disk, are already aligned to that of the 
binary. In this case any information on the initial inclination 
would be lost and the system would evolve as a coplanar case 
dMarzari & SchollL l2000t iThebault etail I2006L [2008L 120091: 
Paardekooper et al.Ll2008h . 

Apart from the case of fast relaxation, which possibly occurs 
in a minority of cases with very high viscosity, in all other cases 
we expect no or very weak coupling between the gas disk and 
the planetesimals orbital evolution for binaries with inclination 
i > 10°. Note th at this low-; case w ith gas drag has been recently 
investigated by IXie & Zhou! d2009l) . who showed that small in- 
clinations between the binary and a circumprimary disc might 
favor planetesimal accretion as compared to the fully coplanar 
case. 

For our numerical exploration, we will thus make the sim- 
plifying assumption that planetesimals evolve in a gas-free en- 
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vironment: the gas-drag force, which introduces a de-phasing of 
the planetesimal perihelia, does not come into play as in the 2-D 
case introducing a de-phasing of the planetesimal perihelia. The 
evolution of the s warm can be describe d as a pure gravitational 
N-body problem (The bault et all 120061) where the relative im- 
pact velocity steadily increases because of both the de-phasing 
of perihelia and nodes. In the next section we will numerically 
compute the evolution of planetesimal relative velocities. 



3. Numerical procedure 

Planetesimal trajectories are computed within the frame of the 
restricted 3-body problem made of the central star, the compan- 
ion star and massless planetesimals. We use the same code as 
in previous studies of the 2-D case (e.g. iThebault et all l2006i 
l2008ll2009h . since this code is 3-D in essence and can handle 
out of plane perturbers. As already mentioned, the main param- 
eter we are interested in is the evolution of the average impact 
velocities within the population of test planetesimals. To that ef- 
fect, our code has a build-in close encounter search algorithm, 
which tracks down at each timestep all 2-body encounters, al- 
lowing to precisely compute the relative velocity for each colli- 
sion (a "collision" being defined as a close encounter within an 
"inflate d radius" equal to 3 X 10~ 4 AU assigned to each parti- 
cle, see Theb ault et all (120091) for more details). The precision 
we obtain in our relative velocity estimate is of the order of 5m/s 
at 2AU. 

These values of (Av) have then to be interpreted in terms 
of accreting or eroding impacts. The limit between erosion and 
accretion is defined by a threshold velocity v*^^, which de- 
pends on the respective sizes s\ and S2 of the impactors, as well 
as on the value of Q* s \ yS 2, the threshold energy for catastrophic 
fragmentation. Unfortunately, the parameter Q* is very poorly 
constrained and estimates found in the literature can differ by 
up to more than 2 orders of magnitude. We chose here a careful 
approach and consider that Q* is comprised between 2 limiting 
values for a "hard" and "weak" prescription. This will in turn 
result in 2 bracketin g values for v* sI s2 (see the discussion in 
Thebault eTaD,|2006l for more details). 

The initial planetesimal swarm is made of 15000 bodies ini- 
tially set on a 2-dimensional disk inclined by an angle io with 
respect to the binary orbital plane. All the bodies in the disk are 
started on circular orbits with semimajor axis ranging from 0.8 
to 6.5 AU. All the nodal lines are parallel since all the bodies are 
clustered in a disk shape. The mass of the primary and secondary 
stars are fixed to 1M and 0.5M© respectively. The binary's or- 
bital parameters a/,, are chosen as free parameters, a/, ranges 
from 40 to 100 AU; beyond those values the perturbations of the 
companion in the initial phases of planetesimal accretion are too 
weak. The binary eccentricity <?/, assumes different values from 
to 0.4. The inclination z'o varies from 0° to 40°. For larger in- 
clinations the Kozai mechanism strongly inhibits planetesimal 
accumulation, as we will see in the following. 

Our initial model setup is based on the assumption that ini- 
tially the planetesimal swarm form a flat disk, in other words that 
it is a dynamically "quiet" system, with all planetesimals ef ree 
and if ree = 0. For a gas rich environement this choice might be 
justified by the fact that the progenitors of the km-sized planetes- 
imals are coupled to the gas and cannot have large e/ ree fl In the 
present case, this "decoupling" is more delicate to define, since it 



1 although this issue might be more complex than this simple pic- 
ture, depending on how planetesimals are formed (see Discussion in 
Thebault et al.2006) 



Semimajor axis (AU) 

Fig. 2. Circulation period of the nodel longitude as a function 
of semimajor axis for planetesimals started on a disk around the 
primary star. The companion star has a semimajor axis a D = 50 
AU, and eccentricity ei, — and different inclinations between 
the binary orbital plane and the planetesimal disk are consid- 
ered. The values of the circulation period are computed through 
direct numerical integration of the equation of motion (3-body 
problem). 

could either be the consequence of the planetesimals vertical dis- 
persion around a coherent gas disc (case 1), or the consequence 
of the gas disc dispersal (case 2). In each case, the relative timing 
between the arrival of the "initial" kilometre-sized planetesimals 
and the decoupling from the gas is difficult to pinpoint. In a worst 
case scenario, we could have an initial orbital distribution where 
some planetesimals have inclination close to i forced while others 
are still around i - 0. This would introduce a high initial free 
relative velocity component that could not be erased with time 
(contrary to the gas-rich case, see Fig. 10 of Thebault et al.2006). 
This difficult issue clearly exceeds the scope of this paper, but 
our results should probably be taken as a lower estimate in terms 
of inhibition of planetesimal accretion 

4. Results 

4. 1 . Dynamical behaviour, encounter velocities 

When the planetesimals feel the binary gravitational pull, their 
node longitude O starts precessing at a rate which strongly de- 
pends on the individual semimajor axis of the bodies as shown in 
Figf2] At the same time the binary perturbations cause a growth 
of the eccentricity and a de-phasing of the perihelia. In Fig|3] 
we illustrate the orbital distribution of the planetesimal swarm at 
t — 5 x 10 4 and t — 10 5 yr when the binary orbital plane is in- 
clined by 30° respect to the planetesimal disk and the eccentric- 
ity of the binary orbit is 0.2 (the semimajor axis is 50 AU). The 
different timescales of perihelia and node circulation are mani- 
fest in the plot. The different degree of randomization of m and 
Q. makes a difficult task to predict the evolution of the relative 
impact velocity between the planetesimals. 

The planetesimal disk moves to the binary orbital plane 
within one orbital period of the outer planetesimals and grad- 
ually looses coherence as a disk. The nodes are randomized and 
the planetesimal Keplerian orbits take them out of the disk plane. 
In Fig|4]we illustrate this behaviour by plotting the positions of 
the planetesimals at t — 0, when they are still grouped in a disk, 
and at t = 1 x 10 s yr when the randomization has disrupted their 
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Semimajor axis (AU) 



Fig. 3. Distribution of the planetesimal orbital elements at t — 
5 x 10 4 yr (green dots) and t — 10 5 yr (red dots). The binary 
parameters are if, = 30°, ei, = 0.2 and at, = 50 AU. 
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Fig. 5. Average impact rate in the planetesimal swarm as a func- 
tion of time for the binary configuration: a\, = 50 AU,ei = and 
ib = 30". Different curves refer to increasing radial distances 
from the primary star. The randomization of the node longitudes 
lead to a sparse planetesimal configuration leading to a reduction 
of the encounter rate. 
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Fig. 4. 3-D spatial distributions of planetesimals at t = (red 
dots) and t — 1 x 10 5 yr (green dots). The randomization of 
the nodal lines disrupt the coherence of the initial disk starting 
from the outer regions. The disk becomes an extended cloud. 
The binary orbital parameters are ay, = 50 AU, e/, = 0. and 
io = 30°. 



initial spatial configuration. We adopt this timescale since it is 
a hig her limit for th e initial planetesimal accumulation process 
(e.g Xissauertll993h 

The randomization of the node longitudes affects the dynam- 
ics of the planetesimal population in two ways. The first one 
is that the sparser distribution of the bodies in space leads to a 
lower impact rate in spite of the growth in the relative velocity. In 
Fig|5]we plot the impact rate as a function of time and radial dis- 
tance. It shows a gradually declining trend as the nodes become 
more randomly distributed. This trend is more marked at larger 
distances from the primary star where the nodal randomization 
is faster. After 10 5 years the impact rate is reduced by 55% at 1 
AU and it drops down by 94% when r = 6AU, independently 
of the planetesimal size. This percentage can be interpreted as 
the fractional reduction of the impact rate compared to that of 
a planetesimal swarm around a single star. Around t — the 
impact rate of our model is not yet affected by the binary pertur- 



bations and it can be taken as representative of the impact rate 
around a single star. 

The other, and more crucial effect is the progressive increase 
of impact velocities, as can be clearly seen in Fig|6] For val- 
ues of ib < 40°, this increase is due to the combination of the 
particles small free eccentricities e/ r Q and the large inclination 
oscillations induced by the inclined companion star. Indeed, the 
small random horizontal excursion due to e / r brings in contact 
bodies having both i and O values increasingly different with 
time as the node oscillations get tighter. This effect is of course 
more pronounced in the outer regions of the disk, where the pull 
of the companion star is felt more strongly. Note however that 
the steady (Av) increase is observed everywhere in the 0.8-6 AU 
region, it is just the pace of this progressive increase which de- 
pends on radial distance. 

For ib > 40° a fully different behaviour is observed and the 
Kozai oscillations come into play. The eccentricity of the plan- 
etesimal orbits begin to grow while the inclination decreases in 
order to keep the action: 



Vh = ^1 - e 2 pl cos{I pl ) 



(2) 



a constant of motion. The period of oscillation of e 2 pl depends 
on the semimajor axis of the planetesimal orbit and, as a conse- 
quence, large differences build up in the distribution of orbital 
eccentricities of the swarm. This leads to much higher impact 
velocities as is clearly shown in Figj7] 



4.2. Effect on planetesimal accretion 

As illustrated in Fig|6]and discussed above, the increase of im- 
pact velocities is an unavoidable consequence of the node ran- 
domization due to the companion star's perturbation that affects 



2 The existence of a non-zero ep component is unavoidable. Indeed, 
for an unperturbed disc of kilometer-sized planetesimals, equilibrium 
encounter velocities are of the order of the bodies escape velocities, i.e. 
a few m.s~', corresponding to e/ r in the 10~ 5 - 10~ 4 range 
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all the simulated 0.8 to 6 AU region. However, in spite of this 
undesired effect, planetesimals might still undergo accretion and 
form planets. The critical condition is that the timescale for both 
the mutual velocity growth and impact rate reduction should be 
long compared to the accretion timescale. More precisely, (Av) 
have to stay at a low, accretion-friendly value long enough so 
that large planetesimals have enough time to build up. When 
the high velocity regime is reached, the growing objects have 
reached a v* sl s2 value high enough to have accreting impacts de- 
spite higher (Av). Of course, studying this effect in detail would 
require to follow the evolution of the planetesimal size distri- 
bution in addition to their dynamical one. This is unfortunately 
beyond the reach of deterministic N-body codes Q. We shall thus 
adopt here a simplified criterion and consider the time at which 
an averaged (v*) is reached for two cases: a "small planetesi- 
mals" case with s m ,„ = 1km and s max = 10km and a "large plan- 
etesimals" case with s m ,„ = 10km and s mnx = 50km, assuming 
that planetesimal sizes follow a centered Gaussian distribution 
between s,„,„ and s max . These two limiting (v*) values are indi- 
cated by light and dark grey areas in FigO the width of these 
areas being due to our careful definition of v* s i, S 2 as being com- 
prised between two extreme values for hard and weak material 
(see Section 3). 

As can be seen in Fig|6] for our example case with a/, = 50 
AU, e/, — and j'j = 30°, the whole system remains accretion- 
friendly for a population of large > 10km bodies for the dura- 
tion of the simulation, i.e. 10 5 years, a conservative timescale 
for runaway growth. For kilometer-sized planetesimals, the sit- 
uation is less favorable and the a > 5 AU region becomes hos- 
tile to kilometre-sized planetesimal accretion after ~ 5 x 10 4 
years. In these regions, planet growth can proceed only if in a 
few 10 4 years planetesimals can grow large enough to have ac- 
creting impacts in a (Av) ~ 50 - lOOm.s -1 environment. Note 
however that even if planetesimal accretion is possible, it can 
probably not lead to t he same runaway growth as around an un- 
perturbed single star dKortenkamp et all 1200 ll) . Indeed, the in- 
crease of the impact velocity, even if it cannot stop accretion, 
significantly slows it down by decreasing the value of the grav- 
itational focusi ng factor onto grow ing objects, (see the detailed 
discussion in iThebault et al.L 120061) . Past this initial planetesi- 
mal growth stage, it is difficult to predict the evolution of the 
swarm at farther stages when big planetesimals and planetary 
embryos will collide at these high relative velocities. Large rela- 
tive inclination s mig ht be produced within the planetary system. 
Ouin tana et al.l d2002l) have simulated the formation of terres- 
trial planets in a Centauri and found that planets may indeed 
be formed with large mutual inclinations if the orbital plane of 
the binary is inclined respect to that of the planetary embryos. 
However, their simulations start from a coherent and flat disk 
of protoplanets, while the growing protoplanets might have al- 
ready developed significant inclinations by the time they collide, 
according to our scenario. This suggests that the final phase of 
planet formation in inclined binaries may need additional inves- 
tigation. 

The situation is very different in the it > 40° cases where 
Kozai oscillations dominate the planetesimal dynamics. As ex- 
pected, these cases are much more hostile to planetesimal ac- 
cretion. As can be seen in Fig|7](for a/, = 50 AU, e\, — and 
ib = 40°), after 10 5 years the impact velocities reach values be- 



3 The size evolution of a planetesimal population, under the effects 
of accreting, cratering and fragmenting impacts, can only be treated in 
statistical particle-in-a-box codes for which the dynamical modelling is 
necessarily very simplified 
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Fig. 6. Average impact rate in the planetesimal swarm as a func- 
tion of time at different radial distances from the primary star. 
The grey bands show the erosion limit for planetesimals 1-10 
km in size (lower band) and 10-50 km (upper band). The initial 
inclination between the planetesimal disk and the binary orbit is 
h = 30°. 
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h = 40°. 



yond the erosion limit for both "small" and "large" planetesi- 
mals in the whole a > 2 AU region. Like in the non Kozai case, 
the rate at which (Av) grows strongly depends on the radial dis- 
tance. As an example, the > 5 AU region becomes hostile to the 
accretion of 1-10 km bodies after only a few 10 3 years, whereas 
it t akes almost 10 5 years f or this to be true at ~ 1 AU. 

Theb ault et al.l (120061) derived analytically the timescale be- 
fore the onset of large impact velocities between planetesimals 
in planar eccentric binary systems as a function of the binary pa- 
rameters. In a gas free environment they computed the degree of 
perihelia randomization required to give high relative velocities 
and the time needed to reach it. Even if the inclined case may ap- 
pear similar because accretion occurs in a gas free environment 
and the relative velocities grow because of the node randomiza- 
tion, it is not possible to derive a similar analytical expression. It 
is a complex task to estimate how the relative encounter veloc- 
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ity depends on the degree of node dispersion since this is a full 
3-D problem. In addition, in this scenario both nodes and peri- 
helia, once dispersed, contribute to the grow of the planetesimal 
relative velocities. 

5. Role of the binary configuration 

In the previous section we have analyzed in details the dynamical 
and accretional behaviour of a planetesimal population for two 
specific test binary configurations. For pedagogical purposes, we 
chose cases with eb — to focus on the effect of the binary incli- 
nation and, in particular, on the transition to a Kozai dominated 
regime when i/, > 40°. We now explore how these results depend 
on all the binary's orbital parameters ab, eb and i/, (the mass ra- 
tio between the 2 stars being kept constant and equal to 0.5). For 
the sake of clarity, the accretion/erosion scenario for each binary 
configuration will be summarized by a single simplified param- 
eter, the semimajor axis a\ within which accretion is possible for 
the l-10km planetesimal population at a threshold timescale of 
tf ~ 5 X 10 4 yr for the inner zone ranging from 1-3 AU and 
tf ~ 1 x 10 5 yr for the outer zone extending from 4-6 AU. 
We assume that planet formation is possible if (Av) is lower than 
v*si ;S 2 for any t < Atf. Under this condition larger planetesimals 
can form and resist to higher velocity impacts. As an example, 
by inspecting Fig|6]we can say that planetesimals beyond 5 AU 
reach the erosion regime before 1 x 10 5 yr while for 4 AU (Av) 
is still below the erosion limit. The inner region within 3 AU has 
always impact velocities lower than v* J i jJ 2 when t < 5 x 10 4 yr. 

The choice of Atf is somewhat difficult and arbitrary. 
Runaway growth in a planetesimal swarm around a single star 
is supposed to start after about 10 4 yrs while after 10 5 yrs ap- 
proximately 33% of the disk mass is supposed to be in runaway 
bodies, according to Wetherill and Stewart! dl993l) . These values 
cannot be directly applied to the binary case since 1) the bi- 
nary perturbations increase the relative velocities between the 
planetesimals accelerating the erosion rate 2) the collisional fre- 
quency decreases with time because of nodal dispersion. In this 
scenario it is difficult to derive a reliable value for the Atf with- 
out knowing the details of the planetesimal accretion process. 
This would be possible only with a statistica l code like the planet 
building code ( Weidenschilli ng et al.l d 1997b ) which, on the other 
hand, cannot model the perturbations of a companion star. Here 
we cautiously consider a value of Atf for the inner and outer re- 
gion of the planetesimal disk which is larger than the runaway 
growth timescale and should somehow be a good estimate for the 
time required by planetesimals to grow large enough to sustain 
further accretion into protoplanets. 

The outcome of this analysis is shown in the form of 2- 
dimensional maps. In Figs[8]we plot a\ vs. (a&,2&) for eb = 0.0, 
eb = 0.2 and e/, = 0.4. The outcome for ib = 0° is only a ref- 
erence value since for low inclinations gas drag comes into play 
and it must be included in the numerica l model (for these l ow 
inclination case, see the recent study by IXie & Zhoull2009h . It 
is noteworthy that in the analysis of the data it never happens 
that accretion is possible beyond 4 AU and is prevented within 3 
AU. 

When the companion star is on a circular orbit (Fig|8l,), the 
effect of inclination is noteworthy only for small values of a/,. 
The secular period of the nodes are short enough to perturb the 
disk on a timescale comparable to the accretion timescale only 
for ab < 50AU. For these small separations and eb - planet 
formation is always possible in the < 3 AU region, except for 
the Kozai regime at ib > 40° Beyond this point the swarm is 
quickly eroded because of the enhanced relative velocities stirred 



up by the nodal randomization. At small value of cib there is a 
strong dependence of a/ on ib showing that at the origin of the 
shrinking of the planet formation zone there is the randomization 
of the nodes. For binary semimajor axes larger than 50 AU, the 
situation is much more favorable to accretion, which can only 
be stopped, before tf, in the Kozai regime with ib > 40°. This 
means that a/, ~ 50AU is approximately the border value within 
which secular perturbations alone are efficient enough to affect 
planet accretion in the non Kozai regime. 

For more eccentric binary orbits, the randomization of both 
nodes and perihelia combine destructively and push the limit for 
accretion at larger values of a&. As an illustration, FigJS]; shows 
the situation for eb = 0.4. In the black zone the relative veloc- 
ity is beyond the erosion value even at 0.8 AU from the pri- 
mary star (the inner limit for our planetesimal population). In 
these cases the formation of planets, in particular in the outer 
regions of the disc, is strongly inhibited. For inclinations lower 
than ~ 10° some accretion is possible within 1-2 AU. However, 
as already stated, our model for if, < 10° is less accurate since 
gas drag may in this case affect the evolution of planetesimals. 
For ab larger than 70 AU we retrieve the dependence of a; on 
ib and at ab = 90 AU, planetesimal accumulation can only be 
perturbed in the Kozai regime. 

6. Conclusions 

In this paper we explore the effect of high (//, > 10°) binary 
inclination on the planetesimal accretion process. The main out- 
comes of our work are: 

- The gaseous disk and the planetesimals decouple because 
of the forced inclination of the companion star. As a con- 
sequence, planetesimal accumulation should occur in an al- 
most gas free environment for inclinations approximately 
larger than 10°. 

- The progressive randomization of the planetesimal node lon- 
gitudes lead to the dispersion of the planetesimal disk that 
expand into a cloud of bodies surrounding the star. The 
sparser configuration leads to a significant reduction in the 
collisional rate. 

- The onset of large mutual inclinations among planetesimals 
causes an increase of the relative impact velocity that may 
halt the planet formation process. This effect is particularly 
strong for ib > 40° where the Kozai mechanism comes into 
play. Below this value, planetesimal accretion might be pos- 
sible, preferentially in the regions closest to the primary star, 
depending on the value of it 

- The possibility of planetesimal accumulation depends on the 
balance between the timescale of node randomization and 
that of planetesimal accretion. For a binary on circular or- 
bit, the value of ab around which this balance occurs is 
around 50 AU. Within this value the secular perturbations are 
fast enough to induce large relative velocity on a timescale 
shorter than the typical timescale for planetesimal accumula- 
tion. Outside this limit planetesimals have probably enough 
time to growth big enough to sustain high velocity impacts. 

- When the binary is on an eccentric orbit, the randomization 
of nodes and periastra can lead to destructive collisions for 
binary separations up to 70 AU. 

- The dispersion of planetesimals in the nodal longitude, in 
those cases where the accretion is effective and lead to pro- 
toplanets, possibly leads to planetary systems a) that form 
on longer timescales because of the reduction of the accre- 
tion rate b) on highly mutually inclined orbits. 
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Fig. 8. Maps showing the limiting values for accretion a/ as a 
function of [ab, //,). The top plot refers to the case with e& = 0.0, 
the middle plot to <?/, = 0.2, and the lower plot to - 0.4. The 
color coding gives different values of ai (in AU), the limiting 
semimajor axis beyond which planetesimal accretion is possible. 
Each square of the map refers to the lower value of the labels in 
the axes. The cases for //, = 0" do not include gas drag so they 
are only indicative. 



perturbations: the planetesimal disk is coherent during accumu- 
lation into rocky planets and core of giant planets. In this sce- 
nario, planets with a significant inclination respect to the binary 
orbit can form. For intial inclinations larger than ~ 40° the Kozai 
cycles may lead, on a long timescale, a planet into a highly ec- 
centric orbit which, at the same time, is very inclined with re- 
spect to the equator of the primary. This occurs because dur- 
ing the cycle the inclination respect to the binary plane is sig- 
nificantly decreased at the eccentricity peak leading the planet 
far from the equatorial plane of th e primary. This might explain 
the ob s erved orbit of HP 80606b dWu & Murray I (l2003l) : lGillonl 
(l2009h : lPontet~ai1(l2009ln . 
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For binary semimajor axes much larger than those we con- 
sidered in this paper the nodal longitude randomization becomes 
much longer. As a consequence, the first stages of planetary 
formation will probably proceed unaffected by the companion's 



